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ABSTRACT 


The  purpose  of  this  report  is  to  synthesize  technological  concepts  of  frac¬ 
ture  by  making  a  historical  review  of  the  literature  from  1913  up  to  the  present 
time.  The  pertinent  aspects  of  fracture  and  the  development  of  relevant  con¬ 
cepts  of  linear  elastic  fracture  mechanics  derivatives  were  delineated  and 
summarized  for  the  prediction  of  the  critical  length  of  fatigue  cracks.  The 
pertinent  aspects  of  fracture  consisted  of  the  synthesis  of  Ingliss,  Griffith, 
Orowan,  Irwin, and  Westergaard's  relevant  theoretical  concepts.  It  also 
delineates  Boyle’s  analytical  and  experimental  results  of  the  Westergaard- 
Irwin  theoretical  compliance  of  through-the -thickness  centrally  cracked  plate 
and  sheet  for  the  determination  of  plane-strain  (Kjc)  and  plane-stress  (KJ 
fracture  toughness  stress-intensity  parameter  of  high  strength  alloys. 
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SECTION  I 

INTRODUCTION  AND  BACKGROUND 


The  proliferous  performance  requirements  for  aeronautical  and  aerospace 
missile  systems  have  resulted  in  needs  for  higher  strength  structural  alloys 
(Reference  1).  The  specified  requirements  of  the  minimum  possible  weight  in 
metallic  structural  components  of  aircrafts  and  missiles'  rocket  motor  casings 
have  led  to  the  escalation  of  the  use  of  very  high  strength  steel,  aluminum,  and 
titanium  alloys  (Schematic  A  of  the  Appendix).  Examples  are  steels  with  yield 
strengths  greater  than  200, 000  PSI  and  both  ferrous  and  nonferrous  alloys 
having  a  strength-to-density  ratio  greater  than  or  approximately  equal  to 
700,  000  PSl/lb  per  cubic  inch  or  inches.  At  these  strength  levels,  many  struc¬ 
tural  engineering  alloys  have  exhibited  brittleness  in  the  presence  of  small 
flaws  under  a  gross  uniaxial  state  of  an  externally  applied  load.  These  flaws 
may  be  presented  in  the  alloy  as  a  result  of  casting  and/or  shaping  and  forming 

processing.  They  could  have  been  induced  during  the  fabrication  of  thin  or 
thick-wall  structural  components  or  pressure  vessels.  In  addition  cracks  or 
flaws  may  have  developed  in  the  alloy  from  obvious  stress  raisers  (stress 
concentrators)  during  in-service  applications  (Reference  2).  As  a  result,  ser¬ 
vice  failures  in  high  strength  alloys  often  occur  by  cataclysmic  fracture  (brittle 
failure)  rather  than  by  gross  yielding  and  distortion.  As  a  result  two  design 
criteria  have  evolved  as  shown  in  Schematic  B  of  the  Appendix.  In  addition, 
these  brittle  failures  happen  at  stress  levels  well  below  the  macro-elastic  de¬ 
sign  stresses,  which  are  usually  based  on  the  engineering  0.  2  percent  offset 
yield  strength  of  the  material  (Reference  3). 

The  complete  elimination  of  minute  surface  defects  and  ultimately  fatigue 
cracks  in  airframe  and  missile  rocket  casings  structural  alloys  is  not  within 
the  practical  framework  of  the  present  state-of-the-art  of  metallurgical 
technology.  Therefore,  the"fail  safe"  design  concept  has  been  developed  to 
provide  adequate  structural  integrity  and  safety  in  the  presence  of  such  condi¬ 
tions  of  metallic  materials  damage  that  is  exposed  to  external  loads.  Ac¬ 
cording  to  this  concept,  multiple  load  paths  or  crack  arrestors  are  designed 
into  the  structure  so  that  a  propagating  crack  cannot  destroy  the  load-carrying 
ability  of  the  entire  structure  (Reference  4).  However,  even  in  commercial 
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aircraft,  it  is  not  feasible  to  design  all  load  carrying  members  on  a  "fail  safe" 
basis.  For  an  example,  the  landing  gear  of  aircraft  are  usually  single  load 
path  structures.  This  is  made  necessary  by  weight  and  space  considerations. 

As  a  consequence,  the  need  for  exceptionally  reliable  aircraft  and  missile 
casing  metallic  structures  has  brought  to  a  critical  state  the  dire  need  for  an 
intrinsic  material  parameter  that  will  measure  the  resistance  to  cracking  of  a 
high  strength  alloy  for  the  prevention  of  catastrophic  failure  under  tensile  (static) 
and  fatigue  (dynamic)  loading  conditions  in  service.  The  "fracture  toughness" 
of  a  high  strength  metal  or  alloy  has  become  of  prime  importance  as  the  frac¬ 
ture  resistance  parameter  that  appears  to  meet  the  required  measurable  cri¬ 
teria.  This  parameter  is  believed  to  be  an  intrinsic  property  of  the  material. 
Therefore,  if  so,  it  should  ideally  measure  the  alloy's  ability  to  dissipate  in¬ 
creasingly  stored  elastic  strain  energy  as  localized  plastic  flow  at  the  tips  of 
defects  or  cracks  under  externally  applied  tensile  loads.  This  ideal  transfer 
of  energy  would  also  appear  to  reduce  the  possibility  of  the  slow  growth  and 
rapid  propagation  of  cracks  at  stress  levels  below  the  yield  strength  of  alloys 
under  external  loading  conditions. 

It  has  been  demonstrated  (Reference  3)  that  the  fracture  toughness  of  high 
strength  alloys  is  one  of  the  major  determinants  of  the  thickness  and  inherent 
stress  level  that  an  alloy  can  be  reliably  used  in  the  design  of  highly  tensile 
stressed  structures  fabricated  of  high  strength  alloys.  The  degree  of  toughness 
of  an  alloy  has  been  among  the  most  important  properties  contributing  to  the 
success  or  failure  of  specific  aircrafts  and  missiles  (Reference  5).  A  know¬ 
ledge  of  the  fracture  toughness,  fatigue  crack  propagation,  corrosion  resis¬ 
tance,  and  the  strength  remaining  at  any  stage  of  development  of  a  crack  is 
important  in  the  application  of  the  fail  safe  concept  to  structural  components. 

The  need  for  this  knowledge  becomes  more  acute  as  vehicle  performance  re¬ 
quirements  increase  at  the  expense  of  decrease  material  performance  and 
reliability.  As  the  exploration  of  space  continues,  the  consequences  of  alloy 
structural  failure,  as  in  manned  spacecraft,  become  continually  more  serious. 

The  purposes  of  this  study  (Parts  I,  n  and  III)  were  to  synthesize  techno¬ 
logical  concepts  of  fracture  by  making  a  historical  review  of  the  literature  and 
to  determine  (Parts  rv  and  V)  experimentally  if  high  strength  alloys  fracture 
mechanics  concepts  could  be  applied  to  the  prediction  of  fatigue  crack's  slow 
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growth  and  the  onset  of  rapid  propagation  under  plane-strain  conditions  at  the 
tip  of  a  crack.  The  primary  aim  was  to  add  to  the  knowledge  of  fatigue  crack 
propagation  technology  in  high  strength  alloys  selection  and  utilization  in  design. 
The  mode  of  failure  under  which  the  fracture  toughness  parameter  of  an  alloy 
is  deduced  from  fracture  mechanics  concepts  was  compared  with  the  failure 
under  fatigue  crack  propagation  (Part  V)  to  establish  whether  the  use  of  the 
fracture  concepts  were  technically  valid  in  predicting  critical  crack  lengths  at 
the  onset  of  rapid  propagation  under  fatigue  loading  conditions.  The  degree  of 
correlation,  quantitatively,  between  the  fracture  mechanics  concepts  and  fatigue 
crack  propagation  predictions  was  pursued  (Part  V)  to  establish  whether  the 
critical  crack  length  at  onset  of  catastrophic  failure  under  fatigue  conditions 
could  be  predicted  from  fracture  toughness  experimentally  obtained  data  and 
fracture  mechanics  mathematical  concepts  deduced  in  Part  IV. 

The  mode  of  crack  propagation  failure  under  uniaxial  tensile  loading  (static) 
to  obtain  fracture  toughness  experimental  data  compared  to  the  mode  of  failure 
under  uniaxial  fatigue  loading  (dynamic)  were  investigated  in  Part  V  relative 
to  the  critical  crack  length  determination  at  onset  of  cataclysmic  fracturing. 

In  the  past,  it  has  been  deduced  by  investigators,  from  a  pure  mechanics 
point  of  view,  that  fracture  toughness  properties  and  concepts  of  high  strength 
alloys  can  be  used  to  predict  critical  crack  lengths  under  fatigue  crack  propa¬ 
gation  conditions  without  rigorous  experimental  verification  rationalized  by  the 
discipline  of  mechanical  metallurgy  (mechanics  plus  physical  metallurgy).  As 
a  consequence,  there  was  a  need  for  a  more  intensive  evaluation  and  analysis 
from  a  mechanical  metallurgy  rationale. 

The  analytical  and  experimental  studies  of  fracture  toughness  Kjc  param¬ 
eter  ability  to  predict  the  critical  length  of  fatigue  cracks  and  the  observations 
of  crack  propagation  velocities  were  conducted  on  7075-T7351  aluminum  alloy, 
as  shown  in  Parts  IV  and  V. 

Cracks  were  experimentally  propagated  in  aluminum  plate  specimens 
under  static  (tensile  strain  rate  of  .05  inches/min.)  and  dynamic  conditions 
(conventional  fatigue  cyclic  rates  of  2000  cycles/min.).  Stress-intensity  fac¬ 
tors  at  the  tips  of  through-the-thickness  fatigue  cracks  and  critical  crack 
lengths  experimental  data  were  obtained  and  analyzed,  as  shown  in  Part  V. 
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SECTION  II 
THE  PROBLEM 

1.  STATEMENT  OF  THE  PROBLEM 

The  purposes  of  this  project  were  to: 

a.  Predict  the  critical  length  of  fatigue  cracks  under  the  stable  fracturing 
mode  and  to  experimentally  substantiate  the  analytical  prediction  (Parts  IV 
and  V). 

b.  Study  the  critical  crack  length  at  instability  under  fatigue  conditions 
in  plate  specimens  under  plane -strain  (depicted  in  Part  V). 

2.  SIGNIFICANCE  OF  THE  STUDY 

a.  A  chronological  and  composite  historical  review  of  the:  (1)  develop¬ 
ment  of  concepts  of  linear  elastic  fracture  mechanics  (Part  I),  (2)  theoretical 
and  analytical  aspects  of  fatigue  (Part  II),  and  (3)  the  thermodynamics  of 
fracture  (Part  III)  was  documented. 

b.  The  quasi-static  plane-strain  (Kjc)  fracture  toughness  of  7075-T7351 
aluminum  alloy  was  determined  experimentally  on  4-,  5-  and  6-inch  wide  plate 
specimens  and  compared  with  data  in  the  literature  (Parts  IV  and  V). 

c.  The  behavior  of  the  plastic  zone  under  quasi-static  and  stable  crack 
propagation  in  fatigue  was  studied  in  Part  V  by: 

(1)  observing  the  size  of  the  plastic  zone  under  quasi-static  and 
fatigue  testing  conditions 

(2)  observing  whether  there  was  stable  extension  under  static  loading 
conditions 

(3)  determining  experimentally  the  stable  and  unstable  crack  growth 
velocities  under  fatigue  conditions 

(4)  determining  experimentally  under  stable  and  unstable  crack  pro¬ 
pagation  whether  the  process  was  isothermal  or  adiabatic  using  liquid  crystals 
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(5)  determining  the  effectiveness  of  variable  stress  concentration  fac¬ 
tors  at  the  tip  of  a  fatigue  crack 

(6)  comparing  the  predicted  critical  crack  length  in  fatigue  conditions 
with  that  predicted  from  quasi-static  fracture  toughness  parameters. 

d.  The  author  has  compiled  a  relatively  complete  bibliography  on  fracture. 

e.  Although  all  references  listed  in  the  bibliography  are  not  used  in  the 
text  of  this  report,  they  are  listed  for  future  use  by  those  individuals  who  are 
working  on  the  details  of  fracture  in  research  or  technology. 
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SECTION  III 

REVIEW  OF  PERTINENT  ASPECTS  OF  FRACTURE 

1.  DEVELOPMENT  OF  RELEVANT  CONCEPTS  OF  LINEAR  ELASTIC 
FRACTURE  MECHANICS 

a.  Continuum  Fracture  Mechanics  Concepts 

(1)  Ingliss -Griffith  Basic  Theoretical  Crack  Model  and  Concepts  for 
an  Elastic  Solid  (Glass) 


cr 


♦Plane-Strain  Conditions  consist  of  tensile  strains  in  two  directions  and  stress 
in  three  directions  with  respect  to  the  X-Y-Z  Cartesian  coordinates. 

♦Plane-Stress  Conditions  consist  of  tensile  strains  in  three  directions  and 
stress  in  two  directions  with  respect  to  the  X-Y-Z  Cartesian  coordinates. 


Figure  1.  Griffith  Model:  Perfectly  Elastic  Solid  (Reference  6) 
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Griffith  used  the  elliptical  crack  model,  because  in  1913  Ingliss  presented  the 
stress  analysis  solution  for  elliptical  cracks  embedded  in  an  elastic  solid. 

The  elliptical  crack  model  details  are  as  follows  (References  7  and  8): 

a 


Accumulation  of  strain  energy 
to  feed  the  crack  extension  da 


Without  the  crack  in  the  elastic  solid,  the  strain  energy  density  is: 


2 


where: 

cr  =  the  stress  =  the  intensity  of  a  force  reaction 
within  the  elastic  body  on  the  application  of  an 
external  uniaxially  applied  load. 

E  =  the  modulus  of  elasticity  of  the  solid  body. 

Ingliss*  solution  gives  the  accumulated  energy  density  at  the  tip  of  the  crack  as: 


8 


E 


(2) 


where  aQ  =1/2  the  original  crack  length  before  the  accumulation  of  strain 

energy  required  to  feed  the  crack  extension  da.  The  corresponding  surface 
energy  associated  with  the  crack  which  was  created  as  a  result  of  the  trans¬ 
formation  of  strain  energy  to  surface  energy  is: 

4a y  (3) 


where  y  =  the  unit  surface  energy 
Physical  Significance  of  y  : 

(1)  Surface  free  energy  =  y  =  (  ^  ^  - )  T,V  =  f(A,V,T) 
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where: 

dF  = 

change  in  free  energy  (ergs)  as  a  function  of 
temperature  (T)  and  volume  (V). 

dA  = 

2 

change  in  free  surface  area  (cm  )  as  a  function 
of  temperature  (T)  and  volume  (V). 

(2)  At  high  temperature  and  for  small  systems,  the  system  will  attempt 
to  reach  minimum  "surface  free  energy"  in  terms  of  minimum  area 

per  unit  volume  ^  ^  ,  a  sphere 

(3)  This  means  that  a  fine  wire  will  tend  to  contract 

(4)  Applied  load  (P)  tends  to  extend  the  wire 

(5)  At  zero  change  in  length  (X),  the  applied  load  (P) 


just  balances  the  "surface  energy"  force.  This  point  of  equilibrium  of  the 
surface  is  the  experimental  method  of  determining  the  "surface  free  energy"  y 
(ergs /cm2). 


2 

Analogy  Between  Surface  Energy  (^ergs/cm  )  and 

p 

Surface  Tension  (r-  dynes /cm): 


r  Ldxj 

I 

_ I 

—  d  X  — | 


Soap  Bubble  Model 


P 


therefore: 


dE  =  yL  dX  =  PdX 


where: 

dE  =  Change  in  total  surface  free  energy 


therefore: 
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where: 

y  =  surface  energy  in  ergs /cm 
-2-  =  the  surface  tension  in  dynes /cm 

The  criterion  for  crack  extension  is  that  the  decrease  in  stored  strain  energy 
of  the  elastic  solid  is  greater  than  the  increase  in  surface  energy.  As  a  conse¬ 
quence,  the  accumulation  of  the  maximum  strain  energy  at  the  tips  of  the  crack 
is  dissipated  as  surface  energy  in  terms  of  creating  free  surfaces  (cracks)  in 
the  elastic  solid. 

The  equilibrium  length  of  the  crack  is: 


therefore: 


2  2 

d  (  4a  y  -  ir  a  a  ) 
da"  E 


=  0 


(4) 


a  equilibrium 


4y 

2/E 

2~ 

v  cr 

2 


2 


2 

TTQ  a 


E 


2Ey 

7 r  a 


(5) 

(6) 

(7) 


where  a  =  the  strength  of  the  plate  or  the  critical  stress  before  fracture 
c 

a  a  — (8) 

cr 

The  fracture  strength  of  a  plate  is  dependent  on  the  size  of  the  crack  and  is 
expressed  as: 


°c 


cr  c  and  aQ  having  approached  a  condition  of  a  critical  are  criterion  for  the 

conditions  of  unstable  growth  of  a  crack  or  the  onset  of  unstable  crack  propa¬ 
gation  from  a  physical  significance  point  of  view.  Unstable  cracking  is 
"catastrophic  propagation"  of  a  crack  in  an  elastically  stressed  body. 
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In  summary,  the  Basis  of  Fracture  Mechanics  is  as  follows: 

a  equilibrium 


(2)  Griffith-Orowan-Irwin  Basic  Theoretical  Concepts  as  Applied 
to  an  Elastic-Plastic  Solid  (Metals) 

From  an  energy  view  according  to  Griffith  Theory  and  Irwin  and  Orowan 
criteria,  fracture  toughness  of  a  material  is  the  component  of  work  irreversi¬ 
bly  absorbed  in  local  plastic  flow  (P)  and  total  surface  energy  ( y)  to  create  a 
unit  surface  area  of  fracture.  For  glass  and  ceramics  the  general  functional 
tensile  stress  (<7)to  propagate  a  crack  (a)  is  as  follows: 


2yE 


7T  <7 


v  7 ra 


as  aQ  — •  aQ  +  da 


J  Ey/a  (9) 

where: 

y  -  Total  surface  energy  of  fracture  surface  a 
E  =  Young's  modulus 

For  metals,  the  general  functional  tensile  stress  (a )  to  propagate  a 
crack  (a)  is  as  follows: 

cr  «  VE(y  +  P)/a  (10) 

where: 

P  =  Irreversible  work  absorbed  in  local  plastic  flow 
that  was  added  to  the  surface  energy  term  in  the 
Griffith  equation  by  Irwin  and  Orowan.  The  other 
terms  are  as  previously  mentioned. 

Griffith  assumed  a  crack  (a)  will  propagate  when  the  strain  energy  released 
(  §E)  per  unit  of  crack  extension  exceeds  the  surface  energy  of  the  newly 
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formed  crack.  Strain  energy  released  (  SE)  under  plane-stress  can  be 
described  for  elliptical  cracks  as  follows  as  determined  by  Ingliss: 


2  2 

*  n  a  a 

t  E 

where: 

8  =  the  decrease  in  the  strain  energy  in  a  plate  per  unit 
thickness  as  a  result  of  crack  length  2a  under  plane  - 
stress.  The  crack  is  through  the  plate  of  unit  thickness. 
a  =  applied  stress  required  to  propagate  a  crack  (a) 
a  =  1/2  original  crack  length  of  unit  thickness 


(ID 


The  surface  energy  and  irreversible  work  of  plastic  flow  (a)  of  a  crack 
of  unit  width  W  and  2a  units  long  can  be  described  as  follows: 


4  a  (y  +  P ) 

where: 

a  =  Total  surface  energy  associated  with  the  crack  or 
potential  energy  of  the  crack/unit  plate  thickness 
and  associated  work  of  plastic  flow. 
y  =  Surface  energy  per  unit  area  or  surface  tension, 
a  =  1/2  the  original  crack  or  flaw  length. 

Griffith  postulated  that  when  small  cracks  are  present  in  materials,  spon¬ 
taneous  fracture  would  occur  when  the  total  energy  of  the  system  or  structure 
was  unchanged  by  small  variations  of  the  crack  length  as  illustrated  in  the 
following  energy  balance  equations: 


d(SE  +  a) 

da 


=  0 


where: 


a  = 


2a  = 

In  other  words, 


The  decrease  in  strain  energy  in  a  system  with  a  crack 
present. 

The  potential  energy  of  the  crack  surface  traction 
forces  (surface  energy)  +  dissipated  energy  of  plastic 
flow.  Orowan  added  plastic  flow  term. 

The  total  length  of  the  crack  of  unit  width. 

Griffith  predicted  sudden-fracture  subcritically  when 


(12) 


11 


AFML-TR-70-202 
Part  I 


the  energy  release  rate  was  equal  to  the  rate  of  energy  absorption  due  to  the 
creation  of  the  new  surface.  Therefore,  the  previously  mentioned  energy 
balance  transforms  to  the  following: 


d  (  8  £  +  a  ) 
da 


d 

da 


/  2  2  . 

( v  a  a  ) 


+  4a ( y +  P ) 


2  ir  cr 


=  4  (  y  +  P) 


(13) 

(14) 


a 


/i  (y  +  P)E 
V  7 ra 


(15) 


For  the  essentially  brittle  materials  which  Griffith  was  considering,  the  energy 
associated  with  producing  the  new  surface  or  increased  crack  length  of  unit 
width  is  in  essence  the  surface  tension.  For  metallic  brittle  materials,  the 
energy  associated  with  resistance  to  catastrophic  cracking  is  surface  tension 
and  plastic  deformation.  The  plastic  work  is  difficult  to  obtain  under  test  con¬ 
ditions.  The  determination  of  the  surface  energy  of  polycrystalline  materials 
is  relatively  unexplored. 


It  appears  that  the  projected  ultimate  practical  applications  would  be  the 
determination  of  the  fracture  toughness  critical  strength  or  critical  flaw  size 
of  metallic  materials  from  handbook  surface  tension  data  and  determination  of 
the  energy  associated  with  plastic  work  from  a  true  stress-strain  diagram 
test.  However,  surface  energies  of  high  strength  alloys  are  not  known  and  the 
depth  of  penetration  of  plastic  strains  would  be  needed. 


Griffith  assumed  that  a  crack  will  propagate  when  the  decreased  stored 

2  2 

strain  energy  term  (  TL Z-  0  )  exceeds  the  surface  energy  (4a y  )  of  the  newly 

formed  crack.  The  sign  is  negative  of  the  decreased  stored  strain  energy 
term,  because  it  is  the  energy  released  in  propagating  the  crack  in  the  ma¬ 
terial.  The  total  surface  energy  a  of  a  crack  of  unit  width  and  2a  units  long 
should  not  be  interpreted  as  y  (the  surface  energy  per  unit  area)  in  Equation 
12.  The  above  theory  predicts  the  strength  reasonably  well  for  bodies  which 
behave  in  a  quasi  cleavage  fashion  such  as  glass  and  ceramics.  However,  the 
fracture  of  metals  is  accompanied  by  the  dissipation  of  energy  P  by  plastic 
flow  adjacent  to  the  fractured  surfaces.  This  energy  absorption  (P)  is  often 
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at  least  1000  times  the  energy  corresponding  to  surface  energy  ( y )  (Reference 
9)  for  metals.  Therefore,  the  stress  parameter crof  the  work  term  required 
to  propagate  a  crack  according  to  the  Irwin-Orowan  modification  of  the  Griffith 
formula  can  be  written  as  follows  for  metals: 


Where: 

Pe  =  an  energy  absorption  factor  at  the  tip  of  the  crack 
where  the  strain  energy  accumulates  to  feed  crack 
extension,  y  +  P  —  Pe  where  y  is  ignored. 

Plastic  flow  is  always  found  at  the  edges  of  the  fracture  surfaces  of  metals. 
Therefore,  the  P  term  is  the  most  important  parameter  for  metals,  but  is 
next  to  impossible  to  measure  and  analytically  evaluate  in  practice.  As  a  con¬ 
sequence,  the  Fracture  Mechanics  Analytical  Criteria  was  evolved  by  Irwin. 

In  summary,  the  reason  for  the  increasing  speed  of  crack  propagation  ac¬ 
cording  to  both  the  Griffith’s  and  Irwin-Orowan's  equations  is  that  once  a  crack 
has  been  initiated  and  the  crack  grows  in  length,  the  stress  required  for  pro¬ 
pagation  continually  decreases.  In  addition,  when  the  environmental  tempera¬ 
ture  is  above  the  transition  temperature  of  the  metal,  P  is  large  and  the 
stress  crrequired  to  make  the  crack  grow  will  also  be  large.  Below  the  tran¬ 
sition  temperature  of  the  metal,  the  metal  is  brittle  and  P  will  be  smaller. 

The  stress  necessary  to  cause  crack  growth,  therefore  will  be  reduced. 

Orowan  and  Irwin  independently  proposed  that  the  term  in  Griffith’s  equa¬ 
tions  be  replaced  by  (y  +  P),  a  term  which  excludes  the  surface  energy  but 
not  the  energy  of  plastic  deformation,  yielding  (Reference  9): 


2(y +  P)E 


7 r  a 


,  mH 

v  it  a 


P)E 


(17) 

(18) 


A  new  viewpoint  was  introduced  by  Irwin  and  Orowan  who  recognized  that  un¬ 
stable  crack  extension  would  develop  regardless  of  plastic  flow,  if  the  plastic 
strains  tended  to  localize  near  the  boundaries  of  the  crack.  As  a  consequence, 
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work  expended  in  plastic  deformation  (P)  replaced  the  surface  tension  (  y)  as 
the  factor  controlling  fracture  toughness. 

This  was  the  beginning  of  "Fracture  Mechanics."  This  trend  of  thought 
states  that  the  strength  of  a  high  strength  metal  or  alloy  is  a  function  of  the 
load  bearing  capacity,  cr  ,  and  the  defect  size,  a.  Metals  and  alloys  are  plas¬ 
tic  materials,  and  there  is  a  restriction  of  plastic  action  at  the  tip  of  a  notch. 
The  energy  expended  in  propagating  a  crack  equals  approximately  (y  +  P).  The 
plastic  flow  at  the  tip  of  the  crack  is  of  engineering  significance  because  P  is 
approximately  lOOOy  for  metals  (Reference  6). 

The  application  of  the  Griffith  model  to  metals  and  alloys  is  complicated 
by  the  following  (Reference  6): 

1.  Metals  and  alloys  undergo  elastic-plastic  response  at  the  tip  of 
a  crack  or  notch. 

2.  Propagation  energy  is  mainly  a  result  of  plastic  flow  and  not 
the  surface  energy. 

3.  If  a  metal  is  fracture  resistant  or  tough  enough  and  (P  +  /  )  high 
enough,  the  metal  or  alloy  will  go  into  gross  yielding  before  a  crack  will  pro¬ 
pagate  in  an  unstable  manner  causing  catastrophic  failure  or  fracture  of  the 
metal  or  alloy. 

However,  high  strength  metals  and  alloys  with  a  strength-to-density  >  700,  000 
inches  have  not  been  developed  that  will  meet  this  structural  criteria  stipulated 
in  3  above. 

(a)  Orowan’s  Basic  Mechanical  Metallurgy  Contribution  Toward  the 
Application  of  the  Griffith  Theory  to  the  Fracture  of  Metals  and  Alloys 

The  commercially  important  high  strength,  structural  aluminum  alloys 
(cr / p  >  700,  000  inches)  are  subject  to  brittle  fracture.  The  low  temperature 
brittleness  can  be  attributed  directly  to  this  fact.  When  fracture  occurs  in 
polycrystalline  high  strength  alloys  by  trans -crystalline  cleavage,  very  little 
energy  is  expended  in  propagating  the  fractures,  so  they  closely  resemble  those 
that  occur  in  glass  or  other  brittle  isotropic  elastic  solids.  Very  little  funda¬ 
mental  research  has  been  done  on  the  mechanics  of  quasi  cleavage  in  poly¬ 
crystalline  metals.  The  most  comprehensive  study  of  crystalline  cleavage  has 
not  been  done  on  metal  crystals,  but  on  crystals  of  an  ionic  salt  (LiF)  by 
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Gilman  (Reference  10).  Much  about  the  mechanics  of  fracture  involving  cleav¬ 
age  was  done  by  Griffith  on  glass.  There  does  not  appear  to  be  any  evidence 
that  anyone  has  observed  true  cleavage  in  a  face-centered  cubic  metal.  Be¬ 
cause  basal  cleavage  is  a  well-known  phenomenon  in  zinc  crystals,  this  does 
not  mean  that  all  hexagonal  close-packed  metals  cleave  on  the  basal  plane  as 
zinc  crystals  cleave.  Magnesium  does  not  cleave  easily  on  the  basal  plane  or 
any  other  plane,  nor  is  there  any  information  in  the  literature  relative  to  cleav¬ 
age  of  cadmium  crystals  (Reference  10).  The  one  important  group  of  metals 
in  which  cleavage  is  most  frequently  observed  is  the  body-centered  cubic 
metals.  However,  the  alkali  metals  are  body-centered  cubic  metals.  They  do 
not  cleave.  The  cleavage  plane  in  the  body-centered  lattice  is  usually  the  (100), 
although  there  are  examples  in  which  it  has  been  indicated  that  cleavage  along 
(110)  is  preferred  (Reference  11).  Cleavage  fracture  in  crystals  and  polycrystals 
is  often  complicated  by  the  fact  that  plastic  deformation  occurs  around  and 
ahead  of  the  cracks  as  observed  by  Irwin  and  Orowan,  independently. 

However,  Orowan  attacked  the  problem  from  a  mechanical  metallurgy 
point  of  view  as  opposed  to  a  pure  continuum  mechanics  approach,  which  as¬ 
sumes  a  perfect  lattice  structure.  Under  the  conditions  of  high  crack-velocity 
and  low  temperature,  plastic  deformation  tends  to  be  supressed,  and  crystal¬ 
line  cleavage  fractures  are  therefore,  equivalent  to  "brittle  fracture"  ma¬ 
terials  such  as  cold  glass.  However,  there  is  one  important  difference.  In 
metals  and  alloys,  cleavage  tends  to  follow  specific  crystallographic  planes  of 
low  indices.  Whereas  in  glass,  this  condition  has  no  meaning  because  the 
material  is  usually  amorphous  rather  than  being  crystalline.  There  is  some 
evidence  that  glass  can  deform  plastically  at  room  temperature  (Reference  10). 
The  amount  of  this  deformation  during  the  period  of  a  brittle  fracture  is  very 
small.  One  may  assume  that  amorphous  materials  deform  elastically  to  the 
instant  of  fracture  (Reference  12).  Therefore,  the  justification  for  the  sys¬ 
tematic  exclusion  of  the  assumption  of  any  plastic  contribution  to  the  fracturing 
process  is  justified  in  the  Griffith  Theory  for  amorphous  materials.  Orowan 
initially  assumed  a  truly  elastic  solid  on  an  atomistic  scale  (one  which  fractures 
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without  plastic  flow  of  any  kind)  placed  in  a  state  of  tensile  stress  by  the  induced 
stress  ( <7^)  and  fracturing  on  the  vertical  plane  A-A  (fracture  plane),  as 

follows: 


A 


cr  - 

N 

c 

c 

c 

c 

c 

^ — WW1  WM — ^ 

A/W 

)-aaaa|aaaa/v^ 
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D 

- cr 

N 

A 
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Figure  2.  Orowan  Model  of  Lattice  Plane  of  a  Crystal  (Reference  10) 


The  two  vertical  rows  of  circles  represent  a  pair  of  lattice  planes  with  an 
interatomic  distance  spacing  of  d  shown  across  the  fracture  plane  A-A  in 
Figure  2.  Orowan  assumed,  if  there  were  no  flaws  in  the  solid  that  fracture 
would  occur  by  breaking  bonds  of  distance  d  between  atoms  which  face  each 
other  across  the  fracture  plane  A-A.  In  the  model,  the  two  vertical  rows  of 
circles  are  atoms,  and  the  bond  between  the  pair  of  atoms  are  indicated  by 
short  horizontal  lines  with  a  distance  of  d.  The  physical  significance  of  this 
model  is  as  follows: 

(1)  The  bonds  can  be  thought  of  as  springs  that  permit  atom  pairs  to  be 
displaced  toward  or  away  from  each  other  depending  upon  whether  the  applied 
external  force  that  induces  cr^.  is  compressive  or  tensile. 

(2)  For  a  small  force  and  small  strains,  the  atomic  displacements  will 
vary  linearly  with  the  applied  stress  (  cr^).  Therefore,  Hook’s  Law  is 
applicable  as  follows: 


€ 


X 

d 


Ee  = 


EX 

d 


(19) 

(20) 
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Where  E  =  Young's  modulus 

€  =  the  strain 

d  =  the  interatomic  distance  in  the  absence  of  stress 

cT  =  the  stress  normal  to  the  fracture  plane  induced  by  a  force 

N 

Orowan's  model  (Figure  3)  for  large  displacements  shows  that  the  re¬ 
lationship  between  the  displacement  (X)  and  the  applied  stress  <JN  is  not 
linear  (Reference  13). 


Figure  3.  Cohesive  Forces  in  a  Solid  as  a  Function  of  Interatomic 
Distance  (Reference  13) 


Therefore,  Hooke's  Law  is  not  applicable.  The  mechanism  of  deformation  and 
fracture  from  the  above  model  can  be  explained  as  follows: 

(1)  When  the  atoms  of  Figure  2  are  forced  close  together,  very  large 
repulsive  forces  are  brought  into  play  causing  the  magnitude  of  the  compressive 
stress  to  rise  more  and  more  rapidly  with  increasing  deformation  as  shown 

in  Figure  3. 

(2)  When  the  atoms  shown  in  the  model  in  Figure  2  are  displaced  in  the 
tensile  direction,  large  positive  displacements  come  into  play  causing  a  de¬ 
viation  from  linearity.  However,  in  this  case  the  restoring  force  decreases  at 
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a  maximum  point  rather  than  grows  in  effectiveness  as  the  displacement  in¬ 
creases.  This  effect  can  also  be  deduced  from  Figure  3,  which  shows  the  curve 
that  relates  the  stress  induced  and  displacement  on  an  atomic  scale. 

(3)  From  Figure  3,  it  can  be  readily  seen  that  when  the  atomic  bonds  are 
lengthened,  the  restoring  stress  rises  to  a  maximum  and  then  falls. 

(4)  The  value  of  the  maximum  stress  (  is  taken  as  the  theoretical 

fracture  stress. 

(5)  At  the  maximum  stress  (  cr^) ,  a  point  of  instability  is  reached 
where  additional  strain  occurs  under  decreasing  stress.  With  the  normal  me¬ 
thod  of  loading  in  tension,  this  condition  leads  to  fracture. 

From  Orowan's  point  of  view  (in  the  study  of  fracture),  the  interesting 

part  of  the  curve  in  Figure  3  lies  in  the  stress  interval  between  0  and  crXTi, 

Nth 

(maximum  restoring  stress-intensity  induced  as  a  result  of  an  external  applied 
force).  Inside  the  region  of  0  to  crNth,  the  fracture  curve  (Figure  3)  may  be 
approximated  as  a  simple  sine  function  of  the  form  (Reference  10): 


(21) 


a, 


N 


Where: 


=  the  applied  induced  stress  as  shown  in  Figure  2 
cr 

Nth  =  the  induced  stress  at  the  instant  of  fracture  of  the  bond 


X  =  the  change  in  mean  atomic  distance  or  displacement 
X/4  =  the  value  of  X  when  the  bond  strength  equal 


The  work  per  unit  area  (W)  of  the  fracture  surface  which  is  expended  in 
creating  the  fracture  surface  on  the  plane  A-A  in  Figure  2  is: 


W 


•  dX 


(22) 


0 


From  Equation  21  which  relates  the  stress  to  the  displacement 
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(X),  one  may  write  for  large  values  of  X: 


dcr. 


N 


dX 


—  ( 
dX  ^ 


CTNth  sm 


2ttX  n  _  2tt  _  27 rX 

X  )  '  X  °Nth C0S  X 


(23) 


For  small  values  of  X,  cos 


27rX 


Therefore: 


dcr, 


N 


dX 


2tt 

X  ^Nth 


(24) 


However,  as  previously  stated  for  small  displacements  of  X: 


<rM  =  E  -3- 
N  d 


dcr, 


dX 


N_  -  JL 

d 


(25) 


Therefore,  equating  Equations  24  and  25  yields: 


2ttct. 


Nth 


_E_ 

d 


X  ... 


(26) 


Solving  for  in  Equation  26  and  substituting  it  into  Equation  22,  yields  the 


work  per  unit  area  as  follows : 


W  =  cr 


Nth  2  7r 


— —  ergs/cm2 


(27) 


cr£  .d 

When  fracture  occurs,  this  element  of  strain  energy  ( — mn  )  was  assumed 


by  Orowan  to  be  transformed  into  the  energy  of  the  surfaces  (2  y )  that  are 
created  by  the  fracture. 

Therefore: 


2  r  ~- 

Nth  0 

(28) 

E 

°Nth  = 

/2jE_ 

V  d 

(29) 
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y  is  the  surface  energy  (surface  tension),  E  is  Young's  modulus,  and  d  is  the 
mean  interatomic  distance  across  the  fracture  plane  (shown  in  Figure  2  at  zero 
stress). 

The  above  equation  (Equation  29)  predicts  extremely  high  theoretical  frac¬ 
ture  strength  for  isotropic  elastic  solids.  Spretnak  rationalizes  that  the  failure 
to  achieve  the  high  *s  because  perfect  periodicity  of  the  atomic  structure 

does  not  exist  because  of  the  existence  of  structural  imperfection  such  as  line, 
point,  and  surface  defects.  Spretnak  has  indicated  further  that  the  onset  of 
plastic  flow  is  achieved  (using  Figure  3  as  a  reference)  before  the  cr^.^  stress 
is  reached,  because  of  the  imperfect  atomic  structure  of  the  lattice  (Reference 
6).  He  stated  that  "travel  up  curve  (Figure  3)  is  interrupted  by  onset  of  plastic 

flow  because  of  presence  of  defects  (point,  line,  and  surface)."  In  many  solids, 

o  11  2  3  2 

d  =  3  A  ,  Young's  modulus  =  10  dynes/cm  ,  y  =  10  dynes/cm  . 

Therefore: 


2 


10  dynes/cm' 


Tensile  strengths  of  this  order  of  magnitude  are  very  rarely  found.  When  ob¬ 
tained,  they  are  found  in  freshly  drawn  glass  fibers  (Reference  14)  and  mica 
sheets  loaded  so  as  not  to  stress  their  edges  (Reference  15).  According  to 
Orowan,  mica  ordinarily  has  a  tensile  strength  of  30,  000  to  40,  000  PSI. 

Orowan  assumed  that  this  low  strength  was  due  to  cracks  located  at  the  edges 
of  the  sheets.  He  further  assumed  that  the  flat  surfaces  of  the  sheets  should 
be  free  of  flaws  since  they  are  cleavage  surfaces  and  should  normally  be  with¬ 
out  defects.  In  testing  this  hypothesis,  Orowan  made  tensile  test  grips  with  a 
width  smaller  than  the  sheet  width.  In  this  manner,  he  loaded  the  mica  sheet 
specimens  so  that  their  edges  were  not  stressed.  The  result  of  this  experi¬ 
ment  was  an  observed  tensile  strength  of  more  than  400,  000  PSI  (Reference  10). 

Griffith's  original  work  on  glass  fibers  presented  strong  evidence  for  his 
theory.  When  freshly  drawn  soda-glass  fibers  were  tested  in  bending,  tensile 
strengths  as  high  as  900,  000  PSI  were  obtained,  which  is  very  close  to  the 

g 

theoretical  strength  of  glass  (  *10  PSI).  Griffith  observed  that  the  strength 

of  glass  is  less  than  1/100  of  its  theoretical  strength  (  cr^^).  This  discrepancy 
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between  the  observed  and  the  theoretical  strength  led  Griffith  to  postulate  that 
the  low  observed  strengths  were  due  to  the  presence  of  small  cracks  or  flaws 
(a)  in  low-strength  glass.  Griffith  assumed  that  the  theoretical  stress  (  cr^^) 
was  obtained  at  the  ends  of  a  crack  (a)  even  though  the  average  stress  (  <rN) 
was  still  far  below  the  theoretical  strength.  This  assumption  was  based  upon 
the  premise  that  the  extremities  of  cracks  have  the  ability  to  act  as  stress 
raisers.  According  to  this  concept,  fracture  occurs  when  the  stress  at  the  tips 
of  the  cracks  exceeds  the  theoretical  stress  (  <xNth).  When  this  occurs,  the 
crack  is  able  to  expand  catastrophically  (References  10,  14,  and  16).  Orowan 
(References  10  and  17)  calculated  the  stress  at  the  ends  of  a  crack  based  upon 
the  Griffith  theory  using  a  flat  plate  with  an  elliptical  hole  of  the  type  shown 
in  Figure  1.  Ingliss  (Reference  18)  computed  the  stress  and  strain  around  an 
elliptical  through-the -thickness  crack  of  the  type  shown  in  Figure  1.  Orowan 
et  al  expressed  the  theoretical  fracture  strength  as  follows: 


Where: 

2  a  =  the  length  of  the  major  axis  of  the  elliptical  through- 
the -thickness  crack 
=  the  average  applied  stress 

p  =  the  radius  of  curvature  at  the  ends  of  the  elliptical  crack 

Orowan  assumed  that  in  an  actual  metal,  the  radius  of  curvature  (p)  at  the 
tip  of  a  crack  is  equal  to  the  mean  distance  (d)  between  atoms  across  the  frac¬ 
ture  plane  A-A  shown  in  Figure  2. 

Therefore: 


From  Equation  29: 
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Therefore: 


a 

Nth 


2a 


N 


Solving  for  the  normal  applied  stress: 


(32) 


(  33) 


Where: 

cr^  =  the  average  applied  stress  at  which  the  crack  will 
increase  in  length  a 

y  =  the  specific  surface  energy 

2a  =  the  through-the-thickness  elliptical  crack  length 
E  =  Young’s  modulus 

The  above  relationship  is  Orowan's  version  of  the  Griffith  criterion  for  brittle 
fracture.  The  numerical  factor  of  Orowan's  modification 


differ  from  Griffith  criteria 


by  approximately  12  percent. 

Again,  it  is  observed  from  the  two  basic  equations  that  as  the  crack 
length  increases  the  stress  (  crN)  to  keep  it  moving  decreases.  This  signifies 
that  once  the  crack  starts  moving  that  it  is  able  to  accelerate  to  high  velocities. 

The  above  criteria  of  Griffith  and  Orowan  applies  specifically  to  a 
through-the-thickness  flat  elliptical  crack  lying  in  a  flat  plate.  Calculations 
for  cracks  of  other  shapes  in  solids  with  different  material  geometries  and 
more  accurate  atomic  considerations  have  been  made  (Reference  10).  The 
general  functional  relations  shown  in  Equation  9  have  been  confirmed.  There¬ 
fore,  they  can  be  considered  of  general  significance  with  regard  to  "brittle 
fracture. " 
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The  above  relationships  and  comparisons  are  applicable  to^  elastic  solids 
(such  as  cold  glass).  Therefore,  the  strain  energy  -  -  a  - — 

that  is  released  as  a  crack  propagates  is  converted  into  the  surface  energy  of 
the  crack  surfaces  and  the  kinetic  energy  of  the  moving  material  at  the  sides  of 
the  crack.  Orowan  considered  the  additional  energy  term  (P),  (Equation  10) 
for  the  propagation  of  crystalline  cleavage  cracks.  This  term  is  the  plastic 
deformation  (P)  occurring  just  ahead  of  the  crack.  The  metal  in  this  region  is 
effectively  in  a  state  of  very  high  uniaxial  tensile  stress  with  the  axis  of  the 
force  inducing  the  stress  located  normal  to  the  plane  of  the  crack.  The  energy 
of  plastic  deformation  also  comes  at  the  expense  of  the  elastic  strain  energy 
that  drives  the  crack  to  create  surface  energy  (y)  in  the  Griffith  criteria.  If 
the  work  to  overcome  plastic  deformation  is  too  large,  the  crack  may  not  ini¬ 
tiate  the  propagation  stage.  This  implies  that  in  those  crystalline  materials 
which  cleave  a  minimum  velocity  must  be  achieved  before  the  crack  will  pro¬ 
pagate  catastrophically.  If  the  crack  propagation  velocity  is  slow,  then  the 
strain  energy  will  be  absorbed  in  the  form  of  plastic  flow  (slip)  predominately. 
This  argument  appears  to  set  forth  the  premise  that  plastic  deformation  (P)  is 
the  most  important  factor  which  can  be  used  to  measure  whether  a  metal  is 
incapable  of  catastrophic  fracturing  or  vice  versa.  The  basic  problem  is  how 
does  plastic  energy  suddenly  decrease  at  onset  of  running  fracture. 

Orowan's  criteria  are  readily  applicable  to  polycrystalline  aggregate.  As 
an  approximation,  one  may  consider  that  the  average  microcrack  has  a  length 
equal  to  a  grain  diameter  (D).  The  previously  defined  analogous  Orowan  cri¬ 
teria  for  the  increase  in  crack  length  (a)  for  an  induced  applied  stress 


may  now  be  modified  by  assuming  that  the  specific  surface  energy  term  (y) 
can  be  replaced  by  a  corresponding  effective  surface  energy  Pg.  This  latter 
quantity  was  used  by  Orowan  to  take  into  account  within  the  fracture  process 
the  true  surface  energy  and  the  energy  of  plastic  deformation  expended  in 
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forcing  the  crack  through  a  polycrystalline  aggregate.  The  Reed-Hill  version 
(Reference  10)  of  the  Griffith  criterion  becomes 


Where: 

P 

e 

D 

°N 

E 


N 


(y+  P)E 
2a 


the  effective  specific  surface  energy 

the  average  grain  diameter 
the  average  applied  stress 

the  Young's  modulus 


(b)  Metallurgical  Significance  of  the  Orowan  Criteria 


(34) 


The  relation  predicts  that  the  strength  of  polycrystalline  metals  which 
fail  by  brittle  fracture  (cleavage)  should  vary  as  the  reciprocal  of  the  square 
root  of  the  average  grain  diameter  (Reference  10).  Empirical  relations  of  this 
nature  have  been  reported  for  both  zinc  and  iron  (References  19  and  20). 

The  physical  significance  of 


a 

N 


/(  X  +  P)E 
V  2a 


is  that  as  the  diameter  of  the  grains  becomes  larger  and  larger,  a  critical 
diameter  is  reached  above  which  the  stress  needed  to  increase  the  microcrack 
length  becomes  smaller  than  the  stress  to  form  or  nucleate  a  crack  inside  a 
crystal.  When  this  occurs,  it  is  conceivable  that  the  first  crack  that  forms 
causes  the  failure  of  the  specimen. 


By  the  time  the  crack  has  grown  to  a  size  equal  to  the  grain  diameter,  it 
will  be  large  enough  to  pass  through  the  boundary  and  continue  to  grow  unless 
stopped  by  second  phase  particles.  Therefore,  above  the  critical  diameter  of 
the  grains,  the  fracture  stress  is  controlled  by  the  stress  required  to  nucleate 
cracks  in  crystals.  Below  the  critical  diameter,  the  fracture  stress  is  con¬ 
trolled  by  the  stress  required  to  propagate  cracks  through  a  polycrystalline 
aggregate.  Since  cleavage  cracks  are  not  nucleated  until  the  yield  point  of  the 
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metal  is  reached,  fracture  should  occur  almost  simultaneously  with  yielding 
at  diameters  larger  than  the  critical  diameter  (Reference  10). 

Gilman  (Reference  21)  has  summarized  the  above  conclusions  in  the  form 
of  a  simple  schematic  diagram  shown  in  Figure  4  (Summary  of  Experimental 
Results  Showing  the  Effect  of  Grain  Size  D  on  the  Yield  and  Fracture  Stresses 
of  Polycrystalline  Metals,  Reference  21). 

Schematic  stress-strain  curves  are  shown  in  Figure  5  corresponding  to  the 
crack  nucleation  and  propagation  limited  ranges  shown  in  Figure  4.  They  are 
depicted  to  further  rationalize  the  nucleation -propagation  limitations  of  poly¬ 
crystalline  metals  (Reference  10).  Figure  5a  of  Figure  5  (Schematic  Tensile 
Stress-Strain  Curves)  gives  the  shape  of  a  curve  that  is  nucleation -limited. 
Figure  5b  illustrates  fracture  that  is  propagation-limited.  There  is  small  but 
definite  plastic  strain  in  the  propagation-limited  case.  Figure  5c  would  be  in¬ 
dicative  of  a  high  (  y  +  P)  term  with  P  »  y  .  The  high  plastic  deformation 
factor  would  make  the  metal  incapable  of  failing  by  cleavage  fracture,  because 
slip  would  occur  more  easily.  Therefore,  the  fracture  mode  would  be  shear 
(Figure  5c)  as  opposed  to  cleavage  (Figures  5a  and  b). 
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E  ° 
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Figure  4.  Summary  of  Experimental  Results  Showing  the  Effect  of 

Grain  Size  D  on  the  Yield  and  Fracture  Stresses  of  Polycrystalline 
Metals.  (After  Gilman,  J.  J. ,  Trans.  AIME,  212,  1958,  p.  783) 
(Reference  21) 
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(c)  Irwin’s  Basic  Mechanics  Contribution  Toward  the  Application 
of  the  Griffith  Theory  to  the  Fracture  of  Metals  and  Alloys 

Irwin  has  written  extensively  on  conditions  giving  rise  to  brittle  fracture 
of  metals  based  on  a  modification  of  the  Griffith  theory  (Reference  9).  The 
analysis  is  based  upon  the  Griffith  criterion  modified  for  energy  absorption  by 
plastic  flow  as  follows  (Reference  9): 


(35) 


a  =  the  critical  fracture  stress  for  a  crack  of  size  2a 

P  =  the  surface  energy  and  the  energy  of  plastic 

deformation 

E  =  the  modulus  of  elasticity 

a  =  1/2  the  crack  length  of  an  elliptical  crack  along  the 

major  axis 

The  analysis  of  the  Irwin  criterion  of  fracture  mechanics  involves  two 
quantities  defined  by  Wells  as  follows  (Reference  9): 

(1)  Strain  Energy  Release  Rate  (G) 

G  is  the  quantity  of  stored  elastic-strain  energy  released  from  a 
cracking  specimen  as  a  result  of  extension  by  a  unit  area  of  the  advancing 
crack  a. 

(2)  Fracture  Toughness  (G  ) 

c 

Gc  is  the  component  of  work  irreversibly  absorbed  in  local 
plastic  flow  and  cleavage  surface  tension  to  create  a  unit  area  of  fracture. 
(Also  referred  to  as  Surface  Energy  Density,  Fracture  Extension  Force.) 

The  Irwin  concept  of  the  "driving  force  per  unit  crack  front"  or  the  strain 
energy  describing  the  physical  significance  of  the  "Strain  Energy  Release  Rate 
(G)"  and  the  "Fracture  Toughness  (Gc)"  for  high  strength  metals  and  alloys  is 
a  more  generally  applicable  and  a  more  readily  useful  concept  than  the 
Griffith  theory.  The  Irwin  concept  does  not  require  the  assumption  of  the  ab¬ 
sence  of  plasticity  and  there  is  no  complication  involving  fixed  grips.  It  shows 
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that  the  "force  tendency  (G)"  exists  for  any  stress  condition  and  that  "G"  can  be 

determined  from  strain  measurements  near  the  crack  tip.  There  are  various 

formulas  for  "G".  However,  since  applications  of  the  basic  "crack  extension 

force  tendency  concept"  to  fracturing  compared  to  the  Griffith  Theory  are  not 

widely  known,  a  brief  description  of  it  is  delineated  subsequently  in  terms  of 

elastic -plastic  contributions.  In  conjunction,  the  critical  value  of  the  crack 

extension  force  G  will  be  defined  as  the  "fracture  toughness"  or  the  "force 
c 

tendency"  (crack  extension  force)  to  resist  the  onset  of  fast  crack  propagation 
in  high  strength  alloys.  The  strain  energy  release  rate  is  physically  defined 
as  follows  (References  22,  23,  and  24).  Assume  one  has  under  consideration  a 
solid  object  acted  upon  by  external  force,  K,  and  the  body  contains  a  growing 
crack,  a..  If,  during  a  time,  St,  the  work  input  from  this  external  force  is 
8wr;  the  change  is  stored  recoverable  strain  energy  SU;  and  the  difference 
causing  the  change  in  stored  recoverable  strain  energy  from  8  to  8u  is  SwQ 

which  consists  primarily  of  a  conversion  of  energy  to  heat  within  the  solid 
object.  Therefore,  one  can  write  a  thermodynamic  energy  balance  as  follows: 


Sw£  =  Su  +  8wq 


(36) 


Sw„  can  be  expressed  in  terms  of  nonrecoverable  and  recoverable  increments 


of  motion  of  F.  since  the  work  input  SwP  is  reduced  to  a  stored  recoverable 

1  £j 

strain  energy  Su  by  an  amount  equal  to  SWq  . 

Therefore: 


Sw  =  I  f.  ( 8&.P^  +  8£.e^ 

E  ;  I  I  I 


(37) 


ments  of  motion  of  F.,  respectively. 

The  thermodynamic  energy  balance  (Equation  36)  has  the  following  physical 
connotation  in  terms  of  the  energy  converted  to  heat  within  the  body: 


W 


Q 


(38) 
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The  elastic  energy  terms  are  £  F  .  S£.e^and  Su.  They  are  defined  as  the 

i  1  1 

"Crack  Extension  Force  Tendency"  or  the  "Strain  Energy  Release  Rate"  as: 

X  F,  Si, 64  -  Su 

G  =  - - - -  (39) 

Sa 


This  is  the  limiting  case  as  8  A,  the  element  of  new  fracture  area  approaches 

zero.  As  deduced  by  Irwin,  the  only  reason  for  this  expression  to  differ  from 

zero  is  the  configurational  change  of  the  object  due  to  extension  of  the  crack  a.. 

The  G  value  is  assumed  to  depend  only  upon  the  applied  forces,  the  configuration, 

and  the  elastic  constants  of  the  material.  The  omission  of  2  F  8^  ^  from  the 

i  1 

expression  for  the  force  tendency  was  a  matter  of  choice  by  Irwin.  However, 
his  argument  was  that  when  one  analyzes  the  longitudinal  yielding  of  a  tensile 
bar  with  no  crack  present,  the  G  terms  in  the  expression: 


SwQ 

Sa 


iFSt^lFS^ 


Sa 


are  of  negligible  importance  compared  to  2  F.  80. 


Pi 


FSfl 


P  a 


Therefore,  the  significant  component  of  the  "force  tendency": 


8w  =  lF,8iiPt+[lFi8tie4  -Su] 
u  i  i 

is 

8wq  =  F  8/  ^  or  F  =  ^Wp .  (40) 

°  SIP£ 

for  the  case  of  the  plastic  "force  tendency". 

One  may  observe  that  on  a  unit  area  basis,  this  plastic  "force  tendency"  is 
just  the  longitudinal  stress. 


Fundamentally,  a  rate  problem  is  the  primary  concern.  Therefore,  the 
time  rate  of  conversion  to  thermal  energy  is  as  follows: 


F  Si  Pfl- 

St 


(41) 
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The  "force  tendency"  for  a  simple  tension  bar  without  a  crack  in  it  is  depicted 
above  and  one  would  want  to  know  the  corresponding  rate  of  yielding.  This 

Sfi. 

yielding  rate  has  to  be  found  experimentally.  The  familiar  plot  of  — 5 -  as 

ot 

a  function  of  F  shows  the  yielding  rate  increases  from  slow  to  fast  over  a 
relative  small  range  of  values  of  the  "force  tendency. "  However,  work  hard¬ 
ening  effects  were  neglected  in  the  plastic  "force  tendency"  analysis.  The 
value  of  F  for  this  range  per  unit  area  is  the  yield  strength  and  is  properly 
regarded  as  an  intrinsic  material  property. 


Following  an  analogous  procedure  as  delineated  above  for  the  longitudinal 
yielding  of  a  tensile  bar  with  no  crack  present,  various  amounts  of  fracture 
extension  force 

X  F  8£.ei-8u 


may  be  applied  to  a  test  specimen  containing  a  crack  and  the  response  mea¬ 
sured  in  terms  of  time  rate  of  fracture  extension.  Generally,  for  high  strength 
alloys  the  extension  rate  change  is  from  slow  to  fast  for  relatively  small 
changes  of  the  crack  extension  force.  The  critical  G  values  (G  =  Gc)  have  a 
significance  relative  to  resistance  to  fracturing  similar  to  that  of  the  yield 
strength  as  a  measure  of  resistance  to  yielding. 


Irwin  has  shown  that  the  crack  extension  force,  G,  is  proportional  to  the 
square  of  the  applied  stress,  a ,  times  the  length  of  the  crack,  a,  as  one  would 
anticipate  from  the  Griffith  analysis  (References  7,  23,  and  24).  Therefore, 
for  a  central  crack  in  a  large  plate  subjected  to  tensile  stress,  <J ,  normal  to 
crack  direction,  the  crack  extension  force  physical  significance  is  as  follows: 

G  .  -  =  -*P-  M2) 

where  2a  is  the  actual  crack  length  and  E  is  Young's  modulus.  This  is  the 
familiar  Griffith  expression  previously  derived  in  this  review.  Using  the 
two-dimensional  situation  (in  Figure  6a  of  Section  III  l.a.  (3)  in  which  stresses 
are  applied  to  a  plate  containing  a  crack  along  the  X  axis  and  such  that  rxy 
along  the  X  axis  is  a  zero  (Reference  25),  one  can  obtain  a  local  stress  in¬ 
terpretation  of  G.  If  the  Cartesian  position  coordinates  in  Figure  6a  of  the 
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stress  function  are  replaced  by  polar  coordinates,  r,  9 ,  using  one  end  of  the 
crack  as  origin  of  coordinates  (Figure  6a),  as  r  approaches  zero,  the  sepa- 
rational  stress,  cry,  approaches 


2r 


1/2 


(43) 


The  Westergaard-Irwin  analysis  treats  the  local  stress  interpretation  of  G  in 
detail  as  summarized  in  Section  III  l.a(3). 


Irwin's  hypothesis  of  the  mechanism  of  the  fracture  process  agrees  with 
the  Westergaard  model  shown  in  Figure  6a.  Cracks  extend  by  a  process  of 
development  and  joining  up  of  new  fracture  origins  near  the  crack  tips.  The 
comparative  extension  possibilities  of  cracks  of  different  length  and  in  different 
fields  of  overall  nominal  stress,  cr,  indicate  that  the  development  of  new  frac¬ 
ture  origins  proceeds  under  the  influence  of  similar  separational  stress  envi¬ 
ronments,  cr  ,  for  two  different  cracks  when  the  G  for  these  cracks  are  the 

y 

same.  In  summary,  the  physical  significance  of  the  Irwin  criteria  means  that, 
in  general,  for  cracks  of  various  lengths  and  locations  in  any  given  solid  ma¬ 
terial  that  the  value  of  G  measures  the  intensity  of  the  crack-tip  stress  field  so 
long  as  the  influence  of  plastic  deformation  accompanying  fracture  extension  is 
limited  to  the  close  neighborhood  of  the  crack. 

As  a  consequence,  Irwin  justifies  the  omission  of  2  F^S-d.  ^from  the 
expression  of  the  total  force  tendency  equation 


8w.  =  I  f.  8£Pi  +  I  f.8£c*  -8u 

Q  .  i  i  .  i  i 

I  I 


and  consider  the  crack  extension  force  tendency  as 

I  F  Siel-  8u 


G  = 


i  1  1 


2 

tt  cr  a 


8a  e 

The  condition  of  fracture  is  accordingly: 

G  =  G_  (Plane  Strain ) 

I  c 


IN 

Glc=  -f  -  <  1 


(44) 
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Where  Kjc  =  stress -intensity  factor  at  the  tip  of  crack  (ajc) 
under  plane-strain  conditions  or  the  "fracture 
toughness"  for  high  strength  metals  and  alloys. 
v  =  Poisson' s  Ratio. 

The  premise  is  that  a  unit  area  of  cleavage-type  fracture  of  a  given  alloy 
at  a  given  temperature  absorbs  the  same  amount  of  energy.  The  quantity  Gjc 

or  KIc  is  considered  a  basic  material  characteristic  (parameter)  and  conse¬ 
quently  a  possible  design  criterion  to  protect  against  brittle  fracture. 

Most  of  the  experimentation  has  been  directed  to  the  determination  of 
Gjc  of  Kjc  for  various  high  strength  alloys.  Test  conditions  for  determining 
Gjc  must  be  such  that  a  brittle  fracture  results.  If  slow  crack  propagation 
precedes  a  running  crack  some  test  method  such  as  compliance,  electrical 
potential  or  resistance,  acoustical  and/or  photography  must  be  employed  to 
discern  the  area  of  crack  initiation  from  slow  crack  propagation.  To  obtain 
values  of  Gjc  experimentally,  it  is  necessary  to  have  expressions  for  G  in 
terms  of  crack  (notch)  dimensions,  geometry  of  the  test  specimen,  modulus  of 
elasticity,  Poisson's  ratio,  and  stress  field.  Fracture  toughness  equations 
have  been  derived  for  infinite,  semi-infinite  and  finite  specimens  for  plane- 
strain  measurements  predominately.  However,  the  center  notch  specimen  can 
be  used  to  measure  the  plane-strain,  and  plane-stress  fracture  toughness  of 
high  strength  alloys  if  the  optimum  thickness  and  the  maximum  finite  width  is 
selected  to  assimilate  the  infinite  plate's  elastic-plastic  stress  distribution. 

The  practical  and  physical  significance  of  the  measure  of  the  "fracture 
toughness  "of  a  high  strength  alloy  is  essentially  a  measure  of  the  ease  and 
speed  with  which  a  crack  propagates  through  the  damaged  material  under  se¬ 
lected  stress  and  environmental  conditions.  Therefore,  the  mose  direct  mea¬ 
surement  methods  are  those  which  measure  the  actual  crack  length  as  a  func¬ 
tion  of  stress  and  time.  In  the  subsequent  Section  HI  l.a.(3),  a  stress  analysis 
criteria  developed  by  Westergaard  and  supplemented  by  Irwin  to  determine  the 
fracture  resistance  of  metals  and  alloys  is  derived.  The  derived  stress- 
intensity  parameters  are  believed  to  be  material  properties  that  describe  the 
toughness  of  high  strength  metals  and  alloys. 
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(3)  Westergaard-Irwin  Theoretical  Stress  Analysis  Concepts  and 
Techniques  for  the  General  Solution  of  the  Plane-Strain  and 
Plane-Stress  Stress -Intensity  Parameters  of  Cracks  in  the 
Griffith-Orowan-Irwin  Model  for  Metals  and  Alloys 

(a)  Westergaard’s  Approach  (Reference  26) 

Westergaard's  analysis  technique  for  the  general  solution  of  problems  of 
cracks  in  elastic  sheets  is  based  on  a  semi-inverse  stress  function  method 
which  involves  the  following  mathematical  model  and  general  theory  of  elasti¬ 
city.  The  problem  considered  is  an  infinite  panel  loaded  in  the  x-y  plane 
(Figure  6a).  An  airy  stress  function,  <£,  must  satisfy  the  equilibrium  and 
compatibility  equations  in  solving  for  the  normal  and  shear  stresses  acting  on 
a  differential  element  (dx  dy)  of  material  in  the  elastic  stress  field  near  the 
crack  tip  (Reference  25). 


For  equilibrium  and  a  state  of  plane  stress,  <j  =  0,  the  following  is 

z 


applicable: 


Scr 

Srxy 

Sx 

8y 

Scr 

St 

.  y  + 

yx 

8  x 
T 


and  the  compatibility  equation  is: 


where 


=  0 


yx 


yx 


V  (crx  +  cry  )  =  0 


(45) 


Sx2  By2 

In  addition,  the  normal  and  shear  stresses  are  given  in  terms  of  cf>  as: 


(46) 


S  Z(p  S 

8y 2  ;  ay  =  s  x2  ’  T*y 


8  2<fr 
SxSy 


(47) 
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Combining  Equations  45,  46,  and  47  results  in  the  biharmonic  equation: 

. v%-7+ 2li?  +17  1481 

By  defining  <f>  as: 

<f>  =  Re  Z  +  ylm  Z  (49) 

where  Z  is  an  analytic  function  of  the  complex  variable 


C  =  x  +  iy 


and 


dZ 


(50) 


results  in  V  ^ <fi  =  0.  Equation  45  then  reduces  to: 

a  =  ReZ  -  ylm  Z ' 
x 

ay  =  ReZ  +  ylm  Z'  (51) 

r  =  -  y  Re  Z ' 

xy  1 

The  stress  function,  Z  (£ ),  which  solves  the  infinite  panel  problem  of 
Figure  6a  is 


z(£)  = 


(52) 


where  c  is  a  uniformly  applied  tension  stress.  Along  the  line  of  expected 
crack  extension  (i.  e.,  along  the  x  axis,  where  9  =  0,  y  =  0),  the  normal 
stress  in  the  y-direction  is: 


In  addition  (Reference  26), 


a  can  be  expressed  as 


(53) 


(54) 
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Therefore,  combining  Equations  53  and  54  results  in 


K 


a 


2  7rr 


2 

x  - 


(55) 


where  x  =  a  +  r.  Since  the  stress-intensity  factor  describes  the  stress  field 
near  the  crack  tip,  the  following  limiting  process  is  taken: 


K 


lim 
r  — 0 


Zjr  (  a  +  r  )2  ~|l/2 

2a  +  r  J 


a  -J  7TQ 


(56) 


Equation  56  is  an  infinite  panel  solution  for  the  crack-tip  stress-intensity 
factor  K.  However,  it  is  of  practical  interest  to  determine  the  stress-intensity 
factor  for  a  centrally  cracked  panel  of  finite  width.  The  approach  is  to  consider 
a  collinear  array  of  cracks  of  length  2a  subjected  to  uniform  tension,  cr ,  and 
spaced  at  equidistant  W  as  shown  in  Figure  6b.  The  appropriate  Westergaard 
stress  function  that  solves  the  given  problem  is  (References  25  and  26): 


z(C> 


<r  sin 


zL 

W 


[< 


sin 


jlL 

W 


)2-( 


sin 


wa  .2  I  i/s 

W  J 


(57) 


Therefore,  along  the  x-axis: 


cr  sin 


7TX 


<ry  =  Re  Z  = 


\ ,  ■  7TX  .  2  .  .  7T0  .  2  1  I/; 

)  -(.in—)  j 


(58) 


W  '  '  '  '  w 

where  x  =  a  +  r.  The  resulting  stress-intensity  factor  is: 


K  =  l1™  cr  sin  no*'}  r 

r  0  w  L 


(  sin 


2  7r  r 

7r(  a  +  r )  x  2 


1/2 


\  2  /  .  7ra  \  2 

)  -(s,n^r) 


W 


(59) 


or 


K  =  CTy?°  [v7  ,on  -^]'  2  <60» 

Therefore,  Equation  60  is  the  stress -intensity  factor  for  a  finite  centrally 
cracked  panel  loaded  in  uniform  tension.  In  comparing  Equations  56  and  60. 
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the  term  tan  ^  j1^2  is  a  stress  field  correction  factor  for  the  finite 

geometry. 

The  previous  stress  analysis  of  the  crack-tip  stress  field  and  stress- 
intensity  factor  for  the  centrally  cracked  fracture  specimen  configuration  is 
based  on  mathematical  elastic  theory  techniques.  Therefore,  certain  cor¬ 
rections  and  limitations  will  be  placed  on  the  solution  for  the  stress-intensity 
factor  due  to  localized  plastic  deformation  at  the  crack  tip.  Of  primary  in¬ 
terest  is  the  size  of  the  plastic  enclave  or  zone  at  the  tip  of  the  crack.  This 
is  readily  determined  by  applying  a  yield  criterion  and  incorporating  the  stress 
field  equations  given  previously  (Reference  26).  The  yield  criterion  which 
gives  results  in  good  agreement  with  test  results  is  the  Henkyvon  Mises  dis- 
tortional  energy  criterion  (Reference  27).  For  the  state  of  plane-stress  or 
biaxial-stress  condition,  the  yielding  theory  can  be  expressed  as: 


(T  O’  +  (X 

x  y  y 


+  3r 


xy 


2 

YS 


(61) 


where  cr  is  the  yield  stress  as  determined  from  an  axial  test.  Substitution 
of  Equation  5  of  Reference  26  into  Equation  61  results  in  the  following  ex¬ 
pression,  which  defines  the  extent  of  the  plastic  zone  in  terms  of  polar  co¬ 
ordinates  r  and  9. 

K2  2  6  [  ,  ^  ,  .  2  9  1 

r  =  27ra2  cos  T  [  1  +  3s,n  ~2~  J  t62) 

YS 

Along  the  line  of  crack  extension,  i.e.,  9  =  0  degrees,  the  width  of  the 
plastic  zone, 

K2 

W  =  r  =  - ^ -  (63) 

2-rra 

YS 

The  width  of  the  plastic  zone  is  considered  as  an  added  effective  crack  length. 
Therefore,  Equation  60,  when  corrected  for  this  additional  crack  length, 
becomes: 


K, 


j-j  [  W  tan  ( 


7T  a 


W 


2  Wo- 


)] 


1/2 


YS 


(64) 
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Based  on  observations  of  experimental  data,  the  ASTM  Committee  on 

Fracture  Testing  of  High-Strength  Materials  (Reference  28)  recommends  that 

centrally  cracked  panels  be  so  designed  that  cr  <0.8  crvc,  where  cr  is  the  net 

n  Y  o  n 

section  failure  stress. 

(b)  The  Westergaard-Irwin  Compliance  Equation  (Reference  29) 

Westergaard  (Reference  25)  has  provided  a  convenient  two-dimensional 
stress  analysis  in  rectangular  coordinates  of  a  very  large  flat  plate  with  tension 
applied  in  the  y  direction  and  a  system  of  cracks  along  the  x  axis  each  of  length 
2a  and  the  center  at  x  =  0,  +L,  +  2L. .  .as  shown  previously  in  Figure  6b. 

The  pertinent  values  are  given  in  terms  of  an  analytic  function,  of  the  com¬ 
plex  variable  z  =  x  +  iy,  in  Equation  57  and  as  follows: 


In  Westergaard's  notation  Z,  Z,  Z'  are  related  by: 


(65) 


and 


Z 


d  Z 
dz  1 


and  Z' 


dZ 

dz 


Z'  =  Z(z)  is  an  analytical  function  of  the  complex 
variable  z  =  x+iy<£  =  ReZ  +  y  ImZ^  (Reference  26) 

According  to  this  method  the  stresses  are  assumed  to  be  functions  of  x  and  y 
only  and  are  expressed  as: 


CTy=  ReZ+ylmZ1 
<Tx  =  ReZ-ylmZ' 
Txy  =  - y  Re  Z ’ 


and  the  displacement,  v,  in  the  y  direction  for  a  plane-strain  situation 


Ev  =  2  ( I  -  t/2)  ImZ  -  {  I  +  t/  )y  Re  Z 


(66) 
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where  vis  Poisson's  ratio  and  E  is  the  modulus  of  elasticity. 


Irwin  has  shown  (Reference  30)  that  the  above  stress  analysis  is  a  good 
approximation  to  the  stress  situation  in  a  centrally  cracked  sheet  specimen 
(with  vertical  and  horizontal  axes  of  symmetry  of  the  specimen  as  y  and  x  axes). 
The  specimen,  with  width  W  =  L  is  considered  to  be  one  unit  of  the  crack 
system  (i.e. ,  the  plate  with  the  system  of  cracks  is  regarded  as  being  made  up 
of  a  large  number  of  centrally  cracked  specimens  placed  edge  to  edge).  Thus 
the  specimen  edges  are  at  x  =  +L/2  and  the  crack  extends  from  x  =  -a  to 
x  =  +a  (Figure  6a). 


Boundary  conditions  require  the  stresses  cr^  and  ir^  to  be  zero  along 

the  borders  of  the  crack  and  the  stresses  cr  and  ir  to  be  zero  along  the  side 

x  xy  ° 

boundaries.  All  of  these  conditions  are  fulfilled  except  the  condition  when 
a  =0  along  the  side  boundary.  Irwin  has  remedied  the  situation  by 

X 

rewriting  cr  as 

X 


cr  =  ReZ  -  ylmZ'  -  cr* 
x  ox 


where  is  a  constant  stress  adjusted  so  that 

CD 

crx  dy  =  0  at  x  =  ±  L/2 

0 

The  introduction  of  the  constant  cr  does  not  make  cr  zero  everywhere 

ox  x 

along  the  side  boundaries,  but  it  does  reduce  o"x  to  a  small  value.  It  is  as¬ 
sumed  that  this  deviation  from  an  exact  solution  is  smaller  than  the  relatively 
small  errors  introduced  in  practical  application  by  the  departure  from  linear 
elasticity  theory  at  finite  strains. 

Irwin  modified  the  equation  giving  displacement  in  the  y  direction  to  make 
it  applicable  to  a  generalized  plane-stress  situation.  The  Westergaard-Irwin 
equation  is 

Ev  =  2  ImZ  -  ( I  +  v )y  ReZ  (67) 


♦Westergaard  stated  that  an  arbitrary  constant  stress  could  be  applied  in  the 
x  direction  without  disturbing  the  validity  of  the  solution. 
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The  Westergaard-Irwin  equation  shows  an  expression  for  compliance  at  the 
specimen  edges.  The  conditions  at  the  edge  of  the  specimen  are  x  =  +L/2 
so  that  z  =  +L/2  +  iy  where  y  will  be  1/2  of  the  specimen  gage  length. 

The  integration 


gives 


] 


Z 


where  C  is  a  complex  constant  of  integration  to  be  evaluated  later  using  the 
known  conditions  at  a  =  0. 

Making  use  of  the  identities: 


sin 


and 


-sin  '  (  -  im  )  =  i  sinh  1  (m) 


then 


sinh 


ReC  +  il  mC 


so  that 
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also  x  =  +L/2 


is  real  so  that 


ReZ  =  Z  at  x  =  ±  L/2 


The  expression  for  displacement  v  is  found  by  substitution  into  Equation 
67.  Following  Boyle  (Reference  31)  the  compliance  can  be  expressed  in  the 
dimensionless  form  Ev/  cr  W  by  multiplying  the  expression  by  the  modulus 
of  elasticity  and  dividing  by  the  specimen  width  w.  Observing  that  L  =  W, 
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(68) 


Im  C  is  evaluated  by  noting  that  at  a  =  0  the  specimen  corresponds  to  an  un¬ 
notched  specimen  and  Ev/  cr  w  is  equal  to  y/W: 
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Where  is  the  mechanical  compliance  of  a  centrally  cracked  finite 

width  specimen. 
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a 

a.  Infinite  Panel  Subjected  to  Uniform  Tension 


<7 


Figure  6.  Centrally  Cracked  Panel  Configuration  with  Infinite  and 
Finite  Width  (Reference  2). 
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2.  BOYLE’S  ANALYTICAL  AND  EXPERIMENTAL  RESULTS  OF  THE 

WESTERGAARD-IRWIN  THEORETICAL  COMPLIANCE  OF  THROUGH- 
THE-THICKNESS-CENTRALLY  CRACK  PLATE  AND  SHEET  FOR  THE 
DETERMINATION  OF  THE  PLANE-STRAIN  (KJc)  AND  PLANE-STRESS 
(K  )  FRACTURE  TOUGHNESS  STRESS  PARAMETER  OF  HIGH 
STRENGTH  ALLOYS  (REFERENCE  2) 

Boyle  (Reference  31)  has  described  a  method  of  measuring  the  mechanical 
compliance  (which  is  the  extension-load  ratio  of  a  centrally  notched  or  cracked 
sheet  specimen)  and  has  demonstrated  how  this  measurement  of  compliance 
can  be  used  to  give  a  reliable  estimate  of  notch  depth  plus  cracks  extending 
from  the  notch  roots.  The  method  is  based  on  the  fact  that  the  compliance  of 
a  centrally  cracked  sheet  is  a  function  of  the  length  of  the  crack  (more  accu¬ 
rately  a  function  of  the  projection  of  the  length  of  the  crack  on  the  transverse 
axis).  As  a  crack  extends,  the  compliance  increases  (stress  now  being  taken 
as  the  gross  stress  calculated  from  the  dimensions  of  the  uncracked  specimen) 
finally  increasing  without  bounds  as  the  crack  reaches  the  edges  of  the  specimen. 
As  tension  is  applied  to  a  centrally  cracked  sheet  specimen,  the  extension  of 
the  specimen  increases  approximately  linearly  with  applied  gross  stress  until 
a  stress  level  is  reached  at  which  crack  extension  occurs.  As  the  crack  ex¬ 
tends,  the  slope  of  the  load-extension  curve  decreases  (load  as  ordinate)  be¬ 
cause  of  the  increased  compliance  of  the  specimen.  When  the  stress  reaches 
an  appropriate  level,  the  specimen  fails  (the  crack  propagates  rapidly  across 
the  remaining  specimen  width).  The  validity  of  the  calculation  of  critical 
crack  length  depends  on  an  assumption  that  the  compliance  value  calculated 
from  the  slope  of  a  line  drawn  from  the  origin  of  the  load-extension  curve  to 
the  point  of  discontinuity  as  shown  in  Figure  7  is  the  same  as  the  compliance 
value  of  a  specimen  that  contains  an  original  crack  the  same  length  as  the  cri¬ 
tical  crack.  Boyle  proved  this  assumption  valid  (Figure  8). 

In  using  the  method  described  by  Boyle  (Reference  32),  data  from  the  com¬ 
pliance  curve  are  used  to  calculate  a  numerical  value  for  the  dimensionless 
compliance  factor: 

Ev 

o-W 
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E  =  the  elastic  modulus  of  the  material, 
v  =  the  total  displacement  between  the  attachment  points  of 

the  compliance  gage  (which  are  two  inches  apart)  divided  by  2. 
cr  =  the  gross  stress. 

W  =  the  width  of  the  specimen. 

The  numerical  value  for  the  compliance  factor  is  selected  on  the  ordinate 
of  the  calibration  curve  shown  in  Figure  8.  This  calibration  curve  was  con¬ 
structed  by  making  load-deflection  curves  on  a  series  of  3  by  12  inch  slotted 
specimens  of  7075-T6  aluminum  alloy  which  were  loaded  only  in  the  elastic 
range.  The  dimensional  units  used  for  the  calibration  curve  permit  the  same 
curve  to  be  used  for  other  materials  and  any  centrally  cracked  specimens  in 
which  the  deformation  is  measured  over  a  length  equal  to  2  W/3  where  W  is 
equal  to  the  specimen  width.  The  corresponding  value  for  can  be  obtained 
directly  from  the  calibration  curve,  and  this  can  be  used  in  the  equation  for 

plane-stress  tests  in  Figure  9  to  determine  G  and/or  K  . 

c  c 

a.  Plane-Stress  Tests  to  Measure  K 

c 

When  sheet  material  is  involved,  fracturing  at  flaws  is  usually  by  mixed 
mode  or  by  plane-stress  mechanisms.  In  other  words,  there  is  considerable 
shear-lip  development  associated  with  the  propagation  of  the  fractures.  In 
testing  of  precracked  sheet  specimens,  the  critical  plane-stress  condition  oc¬ 
curs  at  the  "onset  of  rapid  crack  propagation"  which  usually  results  after  some 
"slow"  crack  development.  For  the  stress  analysis  to  be  applicable,  it  is  nec¬ 
essary  to  know  the  crack  length  at  the  start  of  rapid  crack  propagation  at  max¬ 
imum  load.  The  compliance  gage  discussed  previously  was  developed  for  this 
purpose.  A  typical  compliance-gage  record  for  the  mechanical  type  of  gage  is 
shown  in  Figure  7.  The  slope  of  the  straight  line  that  is  the  initial  part  of  the 
compliance  curve  is  dependent  on  the  initial  crack  length  as  well  as  on  the  di¬ 
mensions  and  elastic  modulus  of  the  specimen.  Ideally,  the  slope  of  the  line 
from  the  point  of  initiation  of  rapid  fracturing  to  the  origin  would  be  dependent 
on  the  length  of  the  crack  at  the  point  of  instability  as  well  as  the  dimensions 
and  elastic  modulus  of  the  specimen.  If  the  specimen  is  unloaded  before  frac¬ 
turing,  the  trace  of  the  unloading  line  will  not  go  through  the  origin  because  of 
the  slight  amount  of  plastic  deformation  at  the  leading  edge  of  the  crack.  How¬ 
ever,  if  the  slope  of  this  line  is  "adjusted"  to  pass  through  the  origin,  it  is  as¬ 
sumed  that  this  compensates  for  the  effect  of  the  plastic  zone  on  the  apparent 
crack  length. 
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b.  Plane-Strain  Tests  to  Measure  K^c 

During  development  of  the  compliance  gage,  Boyle  (Reference  32)  noted  that 
under  certain  conditions  there  was  a  sharp  discontinuity  or  step  in  the  com¬ 
pliance  curve  at  the  upper  end  of  the  initial  straight-line  portion.  This  step 
(Figure  7)  represented  rapid  crack  propagation  with  no  increase  in  load  until 
after  the  rapid  crack  propagation  had  stopped.  Further  examination  indicated 
that  it  was  representative  of  plane -strain  fracturing  (often  referred  to  as  crack 
opening  mode  or  opening  mode  of  crack  extension)  since  it  was  a  square-type 
fracture  not  associated  with  shear-lip  formation.  The  load  at  the  step  in  the 
curve  is  called  the  pop-in  load,  and  the  phenomenon  is  called  meta-instability 
since  it  is  not  necessarily  continuous.  The  equation  for  calculating  the  plane - 
strain  parameters  for  the  center-cracked  specimens  at  the  pop-in  loads  is 
shown  in  Figure  9. 

According  to  the  definition,  the  plane-strain  fracture  toughness  is  associated 
with  unstable  fracturing  which  produces  a  square  fracture  and  a  step  or  pop-in 
load  (Figure  7). 

A  Boyle  room  temperature  compliance  gage  which  was  fabricated  at  the 
AF  Materials  Laboratory  and  used  by  the  author  in  a  previous  program 
(determination  of  the  fracture  toughness  of  several  high  strength  steels)  is 
depicted  in  Figures  10  and  11.  The  compliance  gage  was  supplemented  with 
a  load-acoustical  pickup  system  to  increase  sensitivity  to  plane-strain  load 
pop-in  and  make  it  readily  discernible  for  gross  stress  calculations  for  plane- 
strain  fracture  toughness  material  evaluation. 

c.  Basic  Analytical  and  Experimental  Problems  Associated  with  the 

Determination  of  KT  and  K 
1C  c 

The  basic  analytical  and  experimental  problems  associated  with  the 
determination  of  Kjc  (plane-strain  fracture  toughness  parameter)  and  Kc 
(plane-stress  fracture  toughness  parameter)  are: 

1.  The  theory  is  based  upon  the  concept  of  a  finite  crack  propagating 
in  an  infinite  plate  under  the  influence  of  a  uniformly  applied  tensile  load  and, 
in  turn,  setting-up  within  an  elastically  stressed  body  the  stress-intensity  of 
a  force  reaction  from  the  tip  of  a  crack  outward  to  the  boundaries,  respectively. 
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Under  these  conditions  the  plane-strain  stress-intensity  factor  (Kjc)  has 
been  expounded  by  investigators  to  be  a  material  constant  independent  of  thick¬ 
ness  of  sheet  or  plate.  Whereas,  the  plane-stress  stress-intensity  factor  is 
dependent  on  the  thickness  of  sheet  or  plate  material. 


2.  There  are  basic  experimental  problems  in  the  laboratory  when  one 

attempts  to  simulate  the  infinite  plate  K.  and  K  stress-intensity  conditions 

1C  c 

from  the  tip  of  a  crack  to  the  boundaries  in  a  finite  laboratory  test  specimen. 


3.  The  accuracy  of  laboratory  measurements  of  KT  and  K  depends 

1C  c 

upon: 

a.  the  adequacy  of  acuity  of  the  starting  crack. 

b.  the  remoteness  of  the  points  of  load  application  from  the 
region  of  crack  propagation. 

c.  the  sufficiency  of  specimen  size  to  provide  the  necessary 
elastic  constraint  around  the  plasticity  developed  at  the  tip 
of  the  crack  prior  and  during  noncritical  crack  propagation 
under  plane-strain  and  plane-stress  conditions. 

d.  the  proper  ratio  relationship  between  crack  length  and  finite 
specimen  width  in  order  that  the  nominal  stress  at  the  tips 

of  the  crack  does  not  exceed  prescribed  maximum  theoretical 
bounds  derived  for  the  finite  specimen,  so  as  to  assimilate 
infinite  specimen  plastic-elastic  conditions. 

e.  a  sensitive  method  of  detecting  the  onset  of  slow  crack 
growth  under  plane-strain  conditions  of  stress,  the  termination 
of  slow  crack  growth  under  plane-stress  and  the  simultaneously 
onset  of  catastrophic  fracturing. 
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DEFLECTION  (  in.) 

Figure  7.  Typical  Load- Deflection  Curve  for  Centrally  Cracked  Sheet 
Specimens  Showing  Pop-In  (Reference  32) 


TO 

w 


NOTE: 

THEORETICAL  CURVE  IS  FOR  WESTERGAARD  S 
EQUATION 


Figure  8.  Calibration  Curve  Derived  from  Compliance  Measurements 
on  a  Series  of  Specimens  Having  Various  Slot  Length 
(Reference  32) 
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FOR  plane-strain  tests 


Initial  crack  length  =  —  W  =  2 o 
3  o 


Width 


>  8 


Thickness  - 

IV  K&l-'’) 

v  =  Poisson's  ratio 


c 7  =  Gross  stress  at  pop-in  = 
E  =  Elastic  modulus 


Load  in  pounds 
W  x  Thickness 


EGt 


67TCT 


a  =  One-half  crock  length  at  pop-in ( usually  a0 1 
oQ  May  be  measured  on  the  fracture 
°YS  Y'e,d  strcngth  of  thc  nnaterial  at  0.2  percent  offset 
EGIc 

“  £  Is  often  neglected  if  it  is  very  small  compared  to  a 


YS 


Sira 


YS 


FOR  PLANE-STRESS  TESTS 


Initial  crock  length  =  —  W  =  2a 
3  0 


For 


EG„ 


Width 


Thickness 


greater  than  16  and  less  than  45 


2  7T0 

a  W  ton  — 


c  c 

a  =  Gross  stress  at  maximum  load 

_  Maximum  load  in  pounds 

W  x  Thickness 


a  =  One -half  total  crack  length  at  maximum  load  + 


2 


Note*.  In  one  common  specimen  design  for  sheet  and  plate  to  3/8-inch  thickness,  the  width  is 
three  inches  and  the  length  is  12  inches. 


Figure  9.  Symmetrical  Centrally  Cracked  Tensile  Specimen  for  Sheet 
and  Plate  (Reference  33) 
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Figure  10.  Boyle  Room  Temperature  Compliance  Gage  Load-Extension 

System  Combined  with  a  Load-Acoustical  Measurement  System. 
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Figure  11.  Detail  View  of  Boyle  Compliance  Gage  and  Acoustical  Pickup 
Crystal 
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SECTION  IV 
SUMMARY 


A  chronological  and  composite  historical  review  of  the  development  of  con¬ 
cepts  of  linear  elastic  fracture  mechanics  from  1913  up  to  the  present  time  has 
been  documented  for  a  reference  as  PART  I  of  this  technical  report. 
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2.  ol  =  critical  stress  of  normal  stress  at  the  onset  of  fracture 
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